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Unit commitment state’s the strategic choice to be prepared in order to define 
which of the accessible power plants should be taken into account to supply 
power. It permits utilities to reduce generation price of power. In this paper, 
the unit commitment problem is elucidated by taking N-1-1 contingency as a 
foremost constraint. The standard N-1-1 contingency takes the loss of 
sequential two components in the network having intervening interval for 
network modifications in the middle of two losses. The crucial objective to 
carry out contingency constrictions is to make certain that the operations of 
power system are adequately strong to unexpected losses of the components 
of the network. The optimal scheduling/allocation of the generating units is 
resolved by taking into account the N-1-1 criterion of contingency. 
By considering the N-1-1 criterion of contingency, the problem results to 
give an optimised model which is a linear model of mixed integer form. 
The linear program of mixed integer is a technique of an operational 
assessment in which restriction is imposed on few variables to be integers. 
Primarily benders decomposition was considered but for the improvement of 
results, the algorithm of branch and cut is presented. IEEE 30 bus system is 
taken into consideration and widespread analysis is accomplished to 
associate performance of the system under N-1-1 criterion contingency. 
The computational outcomes determine the value for taking into concern the 
intervening interval for the adjustments of the system with respect to the cost 
and robustness of the system. Later to the above model reliability assessment 
is proposed to calculate the Loss Of Load Expected (LOLE). This model is 
solved using MATLAB/MATPOWER software. 
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1. INTRODUCTION 


The unit commitment implicates the units to be selected which can provide the system with the 
anticipated load over an essential time interval at least cost. Moreover, offers the spinning reserve which is an 
operative reserve of specified margin taking into consideration of the power systems constraints for instance 
least possible time for on/off, up/down ramping limits and capacity of least/extreme generation. Since the 
probability of multiple components in the system of near failures happening simultaneously is comparably 
less, taking into account the interval for intervening for the adjustment of the operator may result in least 
conservative and more decisions of unit commitment which are economical. In recent period, with stable unit 
commitment decisions, optimum flow of power for N-1-1 constrained contingency could be studied by 
interdiction technique [1]. Analysis of interdiction implicates that the subset of the components of the system 
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are identified such that the loss of load of the system would maximize if interdicted during the worst case. 
Outages due to potentially cascading in the course of the contingencies can be avoided by modelling the 
adjustments of the system optimally. Basically contingency constraints make sure that even during the failure 
of a single component the performance of the system will continue ‘as is’. In the case of failure of multiple 
components during the contingency constraints, it the required that the system performance continues at an 
level which is acceptable with shedding of load of some amount and during the recovery interval overloading 
of system is permitted. Due to presence of contingency in power system there might be an unbalance 
amongst supply and demand. 

Contingency due to N-1-1 specifies that either the unit of generation or a line of transmission has a 
loss which is followed by the adjustments of the system. At an advanced interval of time there is an 
additional loss of a line or a generator in the system. Therefore contingency due to N-1-1 has three various 
interval of time namely, base case followed by primary contingency and later secondary contingency. 
During the normal condition, the system functions with no component failure which is referred to as the base 
case. When there is a failure of a single component then the system functions as the primary contingency. 
At t,, there is a change in state from nominal operation to operation due to contingency and primary loss is 
also observed. During the secondary contingency, operation of the system with the failure of two components 
is observed and it’s acceptable for precise shedding of load and overloading of line [2]. In references [3, 4] 
consideration of unit commitment with N-1 constrained contingency is taken into account and employed 
switching of line to ease congestion so as to produce resources of generation which can be dispatched 
economically. In recent times [5], advancement of a robust model of unit commitment has been taken place 
depending on N-k constrained contingency and ambiguity of interval in view of uninterrupted wind/load. 
The solution of the model is obtained by combination of modified method of benders decomposition and 
algorithm of generation by column and constraint (C&CG). As the demand and renewable penetration has 
increased, the importance for analysing the contingencies has also increased. In few circumstances there is 
less redundancy due to the increase in demand. In [6], a new criterion N-k-€ was presented where the demand 
of the system must be met by no less than (1- €/) fraction following the component of system failure of j. 
Two different algorithms were proposed namely benders decomposition and contingency screening 
implicitly. Reference [7], ambiguous sources of renewable energy with a framework of probabilistic failure 
unit commitment constrained is described where for stable set of outages for transmission and generator, 
the probabilities of failure are exhibited as unambiguous functions for the variables of binary scheduling and 
on loss of load expected the penalties are imposed. For day ahead market, unit commitment consisting of 
reliability constrains is presently [8] for combined energy and scheduling of reserve where if there is some 
instantaneous loss in the unit generation then demand needs to be fulfilled with certain probability. 
A stochastic model of chance constrained two stage for solving the problem of unit commitment with two n-k 
criteria for security is presented. Reference [9] presents optimization framework for two stage under 
contingency N-k constraints for combined energy and scheduling of reserve. Benders decomposition 
approach with inequalities for tighter formulation is presented. Previously, Benders decomposition was one 
of the solution approach used to solve the linear programming for mixed integers model for unit 
commitment. The foremost concept of bendersdecomposition is the formation of a main problem and a sub 
divided problem from the problem which sequentially gets divided into sub problems which are independent. 
It solves the linear program of mixed integer in an algorithm which is an iterative process because extensive 
form thru direct solution is not practical. The benders decomposition has a drawback that is execution of the 
algorithm is easy but its computational controllability is restricted since the main problem needs to be solved 
at each iteration so as to select the unit commitment candidate decision. Therefore to overcome this drawback 
an algorithm known as branch and cut is outlined [10]. In [11], one of the approach to solve linear program 
for mixed integers is branch and cut and its principle is, a solution containing smaller subset is obtained from 
a total feasible solution set. Evaluation for these smaller subset can be performed systematically. For larger 
problems of combinatorial the computational time reduces. Unit commitment problem is solved in an hourly 
manner incorporating various constraints of transmission lines and units of generation, uncertainty in 
forecasting of load, random outages of components into reliability problem. In the stochastic unit 
commitment, one of the constraints expected loss of load is included so as to calculate the cost required for 
the reserve to be supplied [12]. The availability of capacity of power at a load bus is assessed using reliability 
calculation in a composite system [13]. 

In this paper, N-1-1 contingency for solving the unit commitment problem is considered. The vital 
difference between this model and previous models is the consideration of the time for intervening in the 
middle of the primary and secondary contingencies. The unit commitment problem with contingency N-1-1 
constraints is solved using an algorithm of branch and cut. Lastly, reliability evaluation is proposed to the 
above model. The LOLE is calculated. Reliability evaluation is carried out for two reasons. Primarily, long- 
term assessments are carried out to support in the planning of the system. Next, short-term valuations help in 
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day to day decisions for operation. The model proposed here for reliability analysis based N-1-1 contingency 
constrained unit commitment has not been analysed previously and this model has advantage over the other 
models because of the intervening time for system adjustment and additionally the adequacy evaluation for 
generating capacity is also evaluated.This manuscript is organised as follows: Initially in this section the 
introduction and literature survey is discussed. Next in Section 2, the mathematical formulation for 
contingency N-1-1 constrained unit commitment is discussed. Section 3 presents the reliability evaluation. 
In Section 4 the results are analysed and in Section 5 the conclusion for work done is discussed [14]. 


2. MATHEMATICAL FORMULATION 
As a part of mathematical formulation, under the circumstances of the contingency scenario, 
the decisions of unit commitment is made. The mathematical formulation is given below [10, 15]: 


2.1. Formulation for Baseline Unit Commitment 


epee CO Fer) (1) 
Such that xeX (2) 
Hp'+ Aft=d', vt (3) 
B,(6{ — Of) — ff = 0, Ve= Gi, j), t (4) 
Ifi< f, vt (5) 
PeXg <5. Pek Veit (6) 
(xq) < pi — pf 1 <ry(x,), Vg,t (7) 


Where, B, : Line e’s susceptance 

f : Capacities of a line vector 

d' : Demand vector in interval t 

Pg »Pg : Lower and upper limits for the capacity of a generator g 


Tg(X),Fg(x): Ramp down/up rates vector for the vector of the committed units x 


H and A: Incidence matrix 
pgfe.OF: The level for generation for g unit, flow of power on line e and the angle on bus 1 in interval t when 
there is no contingency 

In objective function i.e., in (1) the sum of the cost of the start-up and cost for shutdown cS(x) and 
cost for generationcP(p) is being minimized. Here ‘p’ is prescribed as the dispatch levels and for thermal 
generation units the cost for generation is generally specified as quadratic function. Constraint (2) gives the 
information regarding the lowest uptime and downtime for a generators necessity and the cost for start-up 
and shutdown is calculated as the function for the committed units. This term X symbolizes the set of 
possibility. Constraint (3) gives the steadiness of power at bus individually. Constraint (4) gives information 
regarding the flow of power on a line. Constraint (5) gives the limits for the capacity for the line of 
transmission. Constraint (6) denotes the limits for lower and upper dispatch of generator. Constraint (7) 
denotes the limits for ramping of a generator between the two sequential intervals of time. The reformulation 
for the linear program of mixed integer can be obtained by linearizing in piecewise manner the quadratic 
function of cost for generation for the problem of baseline unit commitment. 


2.2. Contingency Modelling 
CafeetO ert (8) 


vinee Ce LeceCe SL, VEE T (9) 
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Leer Cg S LVEEG (10) 
Mrerci S1,Veee (11) 
deter vee c + Liter dete Ge = 2} (12) 


Where, G: Generator set, indexed by g 

E: The transmission line set 

T: The time period set, indexed by t 

Here, C is a set for all contingency N-1-1 situations taking into account all the possible interval of 
periods for the losses in primary and secondary and units of generating and lines of transmission. As c EC, 
following identifications are made, firstly t,;and t, are the two intervals of time denoting the period for the 
losses of primary and secondary respectively. Secondly, ce =1lor ea = 1and ee =1lor ce =1 
represents a component which is failed respectively in individual time period. In the set of the constraints (8)- 
(12), states the failure of two individual components in different period of time (12) where t1 and t2 are the 
initial failure period and failure for second period. Constraint (9) denotes that maximum failure of one 
component can take place in any period of time. Constraints (10) and (11) denotes that the failure of 
individual component happens maximum one time. Constraint (12) requires the failure of two components 
precisely. 

Precise shedding of load and overloading of line are allowable during the period for secondary 
contingency. To ease the infeasibility of operation the factor for overloading of line 0, and allowable 
shedding of load € can be employed. Three different quantities are presented that is ‘wr’ which indicates in 
contingency, ‘h’ denotes the loss of load which is allowable and ‘o’ denotes the overloading of line. 


wt = Yit_ictt, vet 2) 
We = Dierce: Veit (14) 
0 Vit=1,..,tS—1 
hgt = a : (15) 
d; Vit= tt S5an) T 
0 vet=1,..,t$—1 
at -| = (16) 
0, Ve,t=ts5,..., T 


Constraints (15) and (16) needs the loss of load allowable and the factor for overloading are equated 
to djand 0, at the time of secondary contingency. During the remaining interval it has zero value. 


2.3. Contingency N-1-1 Constrained Unit Commitment Model 





min 

x%f,p,8 c°(x) +cP(p) (17) 
f°, p® g°, 0° 
Such that Constraints (2)-(7) (18) 
Hp‘t+ Af + q‘t=d', vt (19) 
B.(0f" — 6f")(1 — we") — fg" = 0, Ve= (i, j),t (20) 
IfSti< f.(1-wS')(1+0S'), Ve,t (21) 
PeXe (1-wg') < pg' S pgxg(1-wg'), Vg.t (22) 
g(%q)(1-wg") S pg — pg? S Ty(Xg)(1-wg"), Vast (23) 
qt < ht, vt (24) 
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Where pests OF are the level of generation for g unit, flow of power on line e, angle on bus ‘i’ in interval t 
when contingency is present. 

The above model is a linear program of mixed integer in tremendously large scale. The objective 
(17) comprises of the generation cost when there is no contingency and the cost of the unit commitment. (19) 
gives the steadiness of power at bus individually with demand q‘°' which is unsatisfied. (20) gives 
information regarding the flow of power on a line. Constraint (21) gives the limits for the capacity for the line 
of transmission with the factor of overloading when contingency takes place at the secondary interval. 
Constraint (22) denotes the limits for lower and upper dispatch of generator. Constraint (23) denotes the 
limits for ramping of a generator. Higher limit on loss of load at bus individually (24). 


2.4. Branch and Cut Algorithm 

Due to the improvement in technologies, the optimization solver made simpler the operation for the 
algorithm of branch and cut. Since the inequalities are added directly within the tree of branch and bound the 
necessity to explore repeatedly the branch and bound tree can be avoided, defined by the variable x of binary 
unit commitment. It includes all the constraints which are dependent of time and doesn’t need any ordering of 
units prior.The relaxation for linear program of main problem be: 


epae Co) + Ch(p) (25) 
Such that xeX (26) 
0<x<l (27) 
(3)-@ (28) 


This technique is used when there is no constraint containing integer value for the solution of linear 
program by means of a simplex algorithm. Suppose a result which is optimum is obtained and instead of 
integer variable there is a value which is non-integer, an algorithm known as cutting plane could be 
considered so as to find the constraints which are linear. At the time of branch and bound, the upper limit of 
linear program relaxation is taken by the solution which is non-integral and the lower limit is served by 
integral solution. In case an upper limit is lesser than present limit which is less, then a node could be pruned. 
The brief algorithm: 

1. Initial integer linear program is added to the active problems list, L 
2. x* is set as null and v* is set as negative infinity 
3. While the list of active problem, L isn’t empty 
a. A problem needs to be selected and removed from L 
Linear program relaxation is solved. 
c. Incase it’s an infeasible solution, head back to 3 (while), or else represent x for solution and v as 
objective value. 
d. Incase v<v’*, head to 3. 
e. Incase if the integer is x, set x* x and v* < v and head to 3. 
f. If desired, then the cutting planes which are disrupted by x are searched. If such any violation, 
need to include them in linear program relaxation and get back to 3b. 
g. The problem is partitioned by branching with feasible regions restriction. Include this problem to 
L, head to 3. 
4. Return x* 


3. RELIABILITY EVALUATION 

Reliability denotes the supply of energy uninterruptedly to end consumers when considering outages 
of scheduled and unscheduled. Due to contingency, the supply and demand might face some unbalance 
between them. Loss of load expected is the anticipation that the supply generation available will not meet the 
load of the system hourly. Also defined as number of days expected that the system would be unsuccessful to 
supply demand of load. LOLE is generally represented in units of time period in place of values which are in 
percentage. If the result for expected loss of load is greater, then it can be concluded that it has deficiencies in 
charging of power or the units which are present are disposed poorly. Calculation for LOLE:The utmost 
appropriate technique is signifying the power system in tabular form which comprises of availabilities, 
un-availabilities and capacities. Generally, a diagram of loads of daily basis is given for consideration. 
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Cumulative probabilities are calculated from the probabilities of combinations of unit outage. Index for 
LOLE is given as a product of the duration of the load expected and the probability of cumulative outage for 
the states essential. LOLE is expressed in hours/year when considering the curve for everyday hourly load 
variation on yearly basis. 


Where p; is the probability of outage for individual capacity and t; is the loss of supply for the duration in 
days. 


4. RESULTS AND ANALYSIS 
An investigation of decisions for dispatch and optimum commitment for IEEE 30 bus system across 
the various criteria of contingency is analysed. In Table | the load profile for 12 hours is mentioned. 


Table 1. Load Profile for 12 Hours 
Hour Demand (MW) 
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Figure 1, 2 and 3 denotes the optimal dispatch levels across various contingencies. The Figure | 
presents the optimum level of dispatch for N-1 contingency. It can be analysed that the generator | is in the 
OFF state during the 1°, 2™4, 3'4, 8, 9! 10", 11" and 12 hour. Generator 2 is in the OFF state during the 1°, 
2nd, gth oth 10", 11 and 12" hour. While the remaining other generators are all in the ON state for all the 
hours and the dispatch levels varies accordingly for each generator. 

The Figure 2 shows the optimum level of dispatch for the N-1-1 contingency. It can be analysed that 
the generator 1 is in the OFF state during the 1%, 2™4, 37? 7%, 8, 9 10", 11" and 12" hour. Generator 2 is in 
the OFF state during the 1‘, 2™4, 8, 9, 10", 11" and 12" hour. While the remaining other generators are all 
in the ON state for all the hours and the dispatch levels varies accordingly for each generator. 
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Figure 3 presents the optimal dispatch level for N-2 contingency. It can be analysed that the 
generator 1 is in the OFF state during the 1“, 2™, 8, 10", 11" and 12" hour. Generator 2 is in the OFF state 


during the 1“, 10", 11" and 12" hour. While the remaining other generators are all in the ON state for all the 
hours and the dispatch levels varies accordingly for each generator. 
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Figure 3. Optimal dispatch level for N-2 contingency 
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Furthermore, the Figures 1, 2 and 3 denote’s the units committed and their dispatch levels in the 12 
hours’ time period. All the units of generator 3, 4, 5, 6 are committed at all hours for the above various 
contingencies. At the 5" hour there is maximum generation and when its 11" hour the generation is the least. 
The dispatch levels varies for different generators at different hour period for different contingencies. 

Next, in Figure 4 comparison is made for total system costs across various criteria of contingency 
for IEEE 30 bus system. The cost for N-1, N-1-1 and N-2 contingencies are $53142, $53361 and $69351. 
It can be concluded that the cost rises uniformly across the contingency constraints N-1, N-1-1 to N-2. 
The difference in the cost across the N-2 and N-1-1 contingency criteria exhibit the value for the time of 
intervening for the adjustment of the system for the duration of a multiple-failure scenario across contingency 


Total system cost $ 
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30000 i — 
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Figure 4. Total system cost under various contingencies 


Lastly, in Figure 5 comparison is made for the expected loss of load across different contingency 
analysis. The loss of load for N-1-1 is 3.788 hr/year and for N-2 it is 3.829 hr/year. As a result, the loss of 
load for N-2 is more compared to N-1-1 because for N-2 simultaneous failure takes place while for N-1-1 
there is an intervening time for system adjustment. If LOLE is more, then the system cost will also increase 
and vice versa. Henceforth from results it can be concluded that for N-2 contingency the loss of load and 
system cost is high as compared to N-1I-1 contingency. In N-1-1, having a time interval for system 
adjustment is better to avoid the losses. 
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Figure 5. Loss of load expected for IEEE 30 bus 


CONCLUSION 
A model of unit commitment taking into account the non-simultaneous failures of components with 


the time for intervening by operators for the adjustment of the system is taken into consideration. This model 
gives a large scale linear program of mixed integer. To solve such model branch and cut algorithm is used. 
The dispatch levels of the system with associated cost of generation across various contingency constraints 
are evaluated and presented. Furthermore to the above evaluation, reliability assessment is presented to 
calculate the expected loss of load for the unit commitment problem on the hourly basis. Comparison is made 
for the results obtained across different contingency and hence it’s concluded that having an intervening time 
for system adjustment is better with respect to cost and expected loss of load. 
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